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Abstract. We discuss the form of the propagator U(t) for the time- 
dependent Schrodinger equation on an asyptotically Euclidean, or, more 
generally, asymptotically conic, manifold with no trapped geodesies. In 
the asymptotically Euclidean case, if x € Cq° , and with T denoting 
Fourier transform, To e~ %T ^ 2t U{t)\ is a Fourier integral operator for 
t 7^ 0. The canonical relation of this operator is a "sojourn relation" as- 
sociated to the long-time geodesic flow. This description of the propaga- 
tor follows from its more precise characterization as a "scattering fibered 
Legendrian," given by the authors in a previous paper and sketched here. 

A corollary is a propagation of singularities theorem that permits a 
complete description of the wavefront set of a solution to the Schrodinger 
equation, restricted to any fixed nonzero time, in terms of the oscillatory 
behavior of its initial data. We discuss two examples which illustrate 
some extremes of this propagation behavior. 



1. Introduction 

Let us consider the Schrodinger initial value problem for a particle moving 
in curved space with metric g 

(1) (D t + ~A + V\il>(z,t) = 0, V(0,-)=V>o, 

with Dt = —idt and A the nonnegative Laplace-Beltrami operator 

— l —d zl g t3 ^gd z j, 

and where V is a smooth, real-valued potential function. The solution is 
given in terms of functional calculus by the formula 

i/>(t, •) = e-^Vo, H = ^A + V, 

and the operator e~ ltH (or its kernel) is called the propagator, or fundamen- 
tal solution. On flat M. n with V = 0, there is an explicit formula for the 
propagator for its Schwartz kernel is 

(2) K(z, w, t) = (2TTity n / 2 e ilz - wl2/2t . 

This may also be thought of as the solution of Q with initial data the 
delta function at w (the initial value problem makes sense for any tempered 
distribution as initial data). 

l 
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This solution exhibits some peculiar properties, from the point of view 
of propagation phenomena. Fixing w = wo € M n , we find that for any 
t > 0, K{t, z, wq) is smooth. Hence the initial delta-function singularity at 
z = wo apparently disappears. Additionally, since the propagator is, by 
self-adjointness of A, unitary on L 2 , we can reverse this process: If we take 

as smooth initial data for a solution to (^Q), then at t = —A -1 the solution 
develops a delta- function singularity, located at z = — £/A. One can see 
this either by explicitly convolving with K(z, w, t) and doing a Gaussian 
integral, or by noting that this initial data is nothing but a multiple of a 
time-translated version of the fundamental solution. 

The upshot, then, is that singularities to solutions of the Schrodinger 
equation can both appear and disappear. In this paper, we address these 
phenomena not just on M n , but more generally on curved space. The basic 
questions are where do singularities go when they disappear, and, conversely, 
what causes their appearance? (These questions are based on the idea that, 
in evolution described by a unitary group, things should not "disappear" 
and "appear," but rather be transformed into something more or less equiv- 
alent.) The most powerful tool for answering propagation questions is to 
understand the kernel of the propagator: a sufficiently precise understand- 
ing of it, particularly its asymptotics at infinity, will allow us to answer both 
questions. We describe a generalization of the formula (j2J and an interpreta- 
tion thereof which enables us to describe the formation and disappearance 
of singularities. We should emphasize that by a singularity at time t we 
mean a point in the wavefront set of ij)(t, •) in the sense of Hormander |Sj. 
We recall that a point in the wavefront set of a distribution is an element of 
the unit cosphere bundle, (z, £) with |£| = 1, which intuitively describes an 
infinitesimal wave located at z with wavefronts normal to £. 

This note is intended as a "user's guide" to the authors' more technical 
paper [7. which describes both a parametrix construction for the Schrodinger 
propagator and, as a corollary, a full propagation of singularities theorem 
answering the two questions posed above. 

Much has long been known about the structure of the propagator on M. n 
with nonzero potentials V; indeed a number of rather precise parametrix 
constructions exist (among others, those of Fujiwara [2], Zelditch ^H]; Treves 
[E|, Yajima [13). 

On curved space, on the other hand, almost nothing was known about 
the propagator until comparatively recently. The first variable-coefficient 
results were those of Kapitanski-Safarov |10| . who showed that in the case of 
a compactly-supported metric perturbation of K n with no trapped geodesies 
the kernel of the propagator is smooth for all t > 0. The same authors 
subsequently constructed a parametrix [2], albeit without the control at 
infinity which will turn out to be essential for the purposes at hand. At 
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around the same time, Craig-Kappeler-Strauss proved the first result about 
microfocal regularity of solutions to (0). Loosely speaking, they showed 
that on asymptotically Euclidean space, if the initial data i[)q is a Schwartz 
function in a conic microlocal neighborhood of a backward geodesic passing 
through a point (z,Q G S*R n , then (z, Q $ WFtp(t) for all t > 0. In 
particular, let x( z ) denote a cutoff function supported in a positive cone 
near infinity. Assume that xWV'o G <S- Now let 7(5) be any geodesic so 
that 7(5) G suppx for all s « 0. We conclude that (7(5), giji^'y (s)) ^ 
WFip(t, z) for any fixed t > 0. Thus a single hypothesis gives regularity along 
the (co-)tangents to a whole pencil of geodesies emanating from suppx for 
all t > 0. In Euclidean space, this yields regularity along the tangents to 
all lines lying in a pencil of directions, and thus constrains the direction in 
which singularities of tp can lie. 

The second author refined the result of Craig-Kappeler-Strauss by in- 
troducing the quadratic wavefront set "at infinity," measuring the quadratic 
oscillation of a distribution 1 . It is analogous to the scattering wavefront set 
previously introduced by Melrose [E] to describe linear oscillations. In 
conditions were given, in terms of the quadratic scattering wavefront set, 
that constrain not only the directions in which singularities can appear, but 
also the times at which they can appear. The actual locations, however, 
remain mysterious from this point of view. Theorems such as these, which 
determine ( and t but not z in the wavefront set (z, Q G WFip(t) have been 
called microglobal results |18j . 

The quadratic-scattering wavefront set results of JH] have recently been 
extended, using a rather different set of tools, to the analytic category by 
Robbiano-Zuily ^3]; in this setting, even defining the quadratic scatter- 
ing wavefront set involves an appropriate two-parameter version of the FBI 
transform, rather than the pseudodifferential methods employed in |121ll6j . 

The first author thanks the Australian Research Council for its support, 
through a Fellowship and a Linkage grant, of this research. The second au- 
thor acknowledges the support of the National Science Foundation, through 
grants DMS-0100501 and DMS-0323021. 

2. The geometry 

The results of [7] hold for a rather general class of manifolds with large 
conic ends, known as scattering manifolds, introduced by Melrose |12| . How- 
ever for simplicity we will restrict ourselves here to a sub-class consisting of 
asymptotically Euclidean spaces. Let z be a Euclidean coordinate and let 

a z 11 
" = TT' r = \ z \ 
\z\ 



^For example, in the case of a distribution e l ^' z - ) , where <f) is homogeneous of degree 
two, hence determined by a function <f> on S n ~ l , the quadratic wavefront set is essentially 
the graph of <j) over the sphere at infinity. 
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be polar coordinates. We will assume that our metric is smooth and has the 
form 

(3) g = ±(dz*f + h, ^^ M^^ fa^o, 

i=i 

where m is a constant and fejj is a C°° function of its arguments, i.e. has an 
asymptotic expansion (Taylor series) in descending powers of r. 

We further make the following crucial nontrapping assumption 2 for g: 

(4) For any geodesic 7(t), lim r(j(t)) = oo. 

t— >±oo 

As for the potential V, we assume V G C°°(R n ;]R) and for r > 0, 

(5) y = £ + ^l 

with c a constant and V a C°° function of its arguments. If m in © and c in 
© both vanish, then H is said to be short-range, otherwise (gravitational) 
long-range. 

We henceforth denote the Hamiltonian 

H= 1 -A + V, 

hence the propagator is e~ ltH . 

3. The form of the propagator 

In order to motivate our results about the Schwartz kernel of the prop- 
agator, let us re-examine the special form © of the Euclidean one. If we 
introduce polar coordinates in the uu variable, we have 

K(t,r,e,z) = e ir2 / 2t (ae-^l*) 

where o = (27rzi) _ ™ //2 e iU ' 2 / 2 *. We have factored K into this form for several 
reasons. First, we are considering K as a function of r, 9 alone, with t > 
and z fixed. Hence a is, from this point of view, constant. Second, while 
K is a rather uninteresting distribution from the point of view of wavefront 
set, it is certainly not a Schwartz function, owing to its oscillatory behavior 
at infinity. We have chosen to exhibit this oscillatory behavior by separating 
out the leading order factor e lv Z 2 *, from the milder, plane wave oscillation 
of e ~ iz 6 '' r / t . Recall now that part of our goal in constructing the propagator 
is to understand the fate of the delta-function singularity launched from 
the point z G 1". From this point of view, the leading order term e tr l 2t 
is useless: it retains no information about where the initial delta function 
lay, and merely records, in its frequency t~ l , the elapsed time. By contrast, 
the plane wave oscillation term is of great interest: we can recover z from 
e -iz-0r/t^ Finally note that t appears in the phase in a very simple way. 

It is actually not necessary to make this assumption, but if we do not then the results 
only apply in the non-trapping part of phase space. 
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Since the quadratic term is of little interest, we will separate it explicitly 
in our formula for the propagator. Thus, for t 7^ 0, and for H short-range, 
let 

W t = e~ ir2 / 2t e~ UH . 

In the Euclidean case, this is a constant multiple of a plane wave e -* 2 ' 6,r /*. 
Fourier transforming in w = r0, then, gives a delta distribution 5(z — w). 
This turns out to be a rather general phenomenon. 

Let T denote the Fourier transform and let x £ C£°(R n ). One of the main 
results from [7j can be expressed as 

Theorem 1. Suppose H is short-range and t ^ 0. Then the operator 

is a classical Fourier integral operator of order on 1" associated to a 
canonical transformation of S*M n . 

Remark. The point of multiplying on the right by x 1S to localize the right 
spatial variable to a compact set. Theorem ^ is essentially concerned with 
the behaviour of the propagator when the left variable approaches infinity. 
What happens as both variables independently approach infinity is more 
complicated. 

Let Q denote the Fourier integral operator of Theorem ^ and 7 the con- 
tact transformation. Then a standard property of FIOs [S] is that Q moves 
wavefront set according to 7: 

WF(Q/) C 7 (WF/). 

It is now easier to see why the propagator creates and destroys wavefront set: 
can can write e~ ttH x = e%r ^ 2t J : ~ 1 Q. Thus Q moves wavefront set around, 
but J-~ l kills it (if it is compactly supported). So a better way to think of 
the above theorem is that Wt maps wavefront set to "Fourier transformed 
wavefront set." 

In the gravitational long-range case, we need to modify the radial variable 
slightly. Let 

m 

(6) r = r + — log r, 

and define 

W t = e- if2 ' 2t e- itH . 

This modification is familiar in, for example, the structure of generalized 
eigenfunctions or of Dollard wave operators for the Schrodinger operator 
with Coulomb potential [2], Jl]. Then To WtX is a quasi-classical 3 Fourier 
integral operator of order zero. 



Quasi- classical here means that the symbol a(x,y,8) has an expansion that includes 
log terms as well as powers of 9. 
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To understand the canonical relation of Theorem ^ will require under- 
standing the long-time limit of geodesic flow. We begin, however, by clari- 
fying the action of Wt itself on singularities of distributions. 

4. Scattering wavefront set 

We have seen that in order to describe the mapping properties of the 
propagator on wavefront set, it is helpful to keep track of Fourier transformed 
wavefront set. The wavefront set of a tempered distribution on W 1 is a closed 
subset of S"*R n ; it is convenient for our purposes to identify this space with 
the product W 1 x 5 n_1 . 

Fourier transforming a compactly supported distribution with nonempty 
wavefront set — a delta-function for instance — yields a smooth function with 
oscillation at infinity — e.g. a plane wave. In tracking where the singular- 
ity has gone under Fourier transform, it is therefore helpful to introduce a 
wavefront set that measures oscillation at infinity. On W 1 , we simply use 
the Fourier transform to do this. 

Definition 2. The scattering wavefront set of a distribution u € 5'(lR n ) is 
the closed subset 

WF sc u cS^'x 

given by 

(2,()eWF scM ^((,i)eWFf( u ). 

So, for instance, 

WF sc e^ z = {(z,0:zeS n - 1 } 
captures the frequency of plane wave oscillation. 

Remark. The scattering wavefront set was introduced by Melrose ^2] i n the 
more global context of scattering manifolds, where the definition is more 
subtle. Melrose's definition encompasses the set defined above, which we 
think of as a subset of an appropriately scaled cotangent bundle at infinity, 
and the ordinary wavefront set inside the cosphere bundle, as well as a third 
component which interpolates between the two, and is a subset, in our non- 
invariant notation, in S™~ 1 x S?" 1 (the cosphere bundle at infinity). 

The quadratic- scattering wavefront set of ^1 can be defined quite simi- 
larly: again a subset of S 11-1 x 1", it can be defined by 

WF qsc u = WF sc u 

where u is the distribution defined by 

u(z) = u(z/^/\z\). 

Theorem ^ now implies that Wt = e - 1T ' 2 l 2t e ~ ltH maps scattering wave- 
front set to ordinary wavefront set, and vice versa. In particular, given the 
non-invariance of the Fourier transform, it is best from a geometrical point 
of view to think of Wt as a "scattering Fourier integral operator" with a 
wavefront relation interchanging W 1 x S 71 ^ 1 and S"™ -1 x R n . 
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5. The canonical relation 

The canonical relation of Wt as a scattering FIO is related to the limit 
of geodesic flow. Given (z, Q G S*R n R n x S 71 ' 1 we let j(t) be the unit 
speed geodesic with 7(0) = z, (7'(0)) 4 = g^Cj, and define (in the short range 
case j 

S f ( z , C) = (0, G S n_1 x R n , with £ = A# + fi, 
7(*) 



(7) 



6 1 = lim . 

t-oo | 7 (i)|' 

A = lim i- 

t— »oo 



m = t lim It(*)I 

t— >oo 



h(t)\J- 



We also define Sft(z,£) = ~Sf{z, —(,)■ These are the forward and backward 
sojourn relations, respectively. We interpret the components of these maps 
as follows: 8 is the asymptotic direction of the geodesic, A is the "sojourn 
time," a measure of how long the geodesic lingers in the finite part of W 1 
before heading off to infinity, and fi measures the angle of contact of the 
geodesic with the sphere at infinity, or, equivalently, distinguishes differ- 
ent geodesies among the pencil of all geodesies with the same asymptotic 
direction. 

We may endow W 1 x S 1 ™ -1 with a contact structure using the canonical 
one- form on S*M. n ; by switching coordinates, as in the definition of scattering 
wavefront set, we may thus endow S"™ -1 x lR n with a contact structure. 

Proposition 3. The maps 

S f ,S b :R n x S"- 1 -> S*™- 1 x R n 

are contact diffeomorphisms. 

Hence Sf and S b are eligible to be quantized to scattering Fourier integral 
operators. 

Theorem 4. The canonical relation of Wt is \t\~ 1 Sf fort < and \t\~ S b 
for t > 0, with the scaling acting in the fiber variable. 

6. A SKETCH OF THE CONSTRUCTION 

In [Jj a result considerably more detailed than Theorem ^ is proven: it 
is shown that Wt is a scattering-fibered Legendrian distribution in both 
space and time variables simultaneously. Legendrian distributions are a 
class of distributions introduced by Melrose-Zworski in scattering theory 
jllj : they are given by oscillatory integrals and on W 1 are essentially the 
same as Fourier transforms of Hdrmander's Lagrangian distributions. The 



4 In the long range case there is a logarithmic divergence in \imt — which needs 

to be removed. We omit the details. 
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scattering-fibered Legendrians are a refinement introduced by Hassell-Vasy 
E] • This finer characterization of Wt amounts to the statement that it is 
a linear combination of distributions of the form 

(8) rt-f ( a(i,r-\M,^e i ^ 1,0 '^ )r/t d?; 

with a and <f> smooth in the short-range case (and having expansions in 
the system r _n (logr) fc for k < n in general), and with cf> satisfying a non- 
degeneracy condition; see [SJ |S] for more details. Note the special role that 
t plays in the phase. All the geometric information is, as usual, encoded in 
the phase function 4>. 

The main step in the proof is the construction of a parametrix for the 
propagator in the class of scattering-fibered Legendrians. One begins this 
process for t near zero near the diagonal of R n x R n (and within the support 
of x). Here an ansatz for the propagator is the WKB expression 

(9) r^/V*^/*^*^ (*,«;). 

i>o 

Solving for the phase in the usual way gives the eikonal equation 

(10) $ = I|V,$| 2 

for $, which has a smooth solution = ^d g (z, w) 2 , where d g is the distance 
with respect to the Riemannian metric g. This ansatz cannot be used for 
z outside the injectivity radius of w, but a more complicated WKB-type 
ansatz 

(11) t -n/2-*/2 f e i *( z ' w > v V t J2t j a j (z,w,v)dv 

J** -> 

remains valid for z in any given compact set. The main problem is to write 
down a suitable ansatz encapsulating the asymptotics as \z\ — > oo. 

In general, one solves a nonlinear first order PDE such as (jlOj) by the 
method of characteristics. In this case the characteristics are given by curves 

s^(z = 7 (s),C< = g i3 {i{s)y,t = 0,r = s 2 /2) 

where 7 is the arc-length parametrized geodesic starting from (w,£) £ S^M™, 
and r is the dual variable to t. The set of all such geodesies sweeps out a 
submanifold L of x M™ x M T which remains smooth beyond the injectivity 
radius. It is a Legendrian submanifold with respect to a naturally defined 
contact structure on R" x x R T . If the projection of L to the z variable 
is a diffeomorphism, then one can write r = &(z) and this determines the 
phase function in @. When this is not true, then one has to use the more 
complicated expression fil[). where one requires k extra variables in order 
to write r as a function of z and v if the null space of the differential of the 
projection to R™ has dimension k. 
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Along the geodesic, the form Q of the metric implies (in the short-range 
case) that 

(12) r( 7 ( s )) = s + sKo + o( s - 1 ), c(7(*)) = rrl + °( s " 1 )- 

|7(«)I 

Hence r( 7 (s)) = r 2 /2 + rZ(w,£) + O(l). Recall that we are interested in 
the phase of the propagator after an e ir / 2< factor has been stripped out. 
Hence defining (p = <E> — r 2 /2, we find that 4> ~ r ^i x iO a t infinity. Since $ 
(or more geometrically, the Legendre submanifold L) has good asymptotics 
at infinity, we can hope to write down an ansatz for the propagator which 
is accurate out to infinity and encodes the oscillations at spatial infinity as 
well as as t — > 0. In fact, this turns out to be possible in the class of fibered- 
scattering Legendre distributions and this gives a parametrix which differs 
from the true propagator by a kernel which is Schwartz in z and 0(t°°), 
together with all its derivatives, as t — » 0. 

Moreover, in the same way that $ in the ansatz is associated to a 
Legendre submanifold, so is 0(0, •, •, •) in the expression (jHJ) for Wt- This is 
proved essentially by symplectic reduction. The Legendre submanifold as- 
sociated to 4> at r = oo then turns out to be the graph of the transformation 
Sf, which is therefore a contact transformation. 

The long-range case is not essentially different. 

In the special case of flat R n , the sojourn time A along the geodesic ema- 
nating from z in direction 9 is 

lim t-\z + t9\= -z-9. 

t^oo 

Note that this is exactly the phase of Wt in the Euclidean case, with no 
oscillatory integral necessary. This is typical of the following special geo- 
metric case: if for each w and 9 ranging over a pair of open sets, there exists 
a unique geodesic 7 with 7(0) = w, lim^oo "y(t)/\'j(t)\ = 6, then we can let 
S(w, 6) denote the sojourn time along this geodesic (i.e. the function S of 
()12|) . but parametrized by different variables). Then S(w,9) parametrizes 
the Legendrian distribution, and we may write the propagator for t > and 
w £ U <s R n simply as 

with a a smooth function. The sojourn time S(w, 9) is closely related to a 
sojourn time defined on pairs of points in S 71 ^ 1 by Guillemin [3] in studying 
the high-frequency asymptotics of the scattering matrix. 

7. Propagation 

As a consequence of the characterization of e ~ ir2 / 2t e ~ ltH as a scattering 
FIO, where f is given by © (f is just equal to r in the short-range case), we 
can now state a propagation theorem, describing precisely when and where 
singularities can appear in R n : 
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Theorem 5. Let ip(t) be a solution to (QJ. Fix a time to. Let (z, () G 5*M n , 
and let 

(9,0 = S b (z,C), (0,i)=S f (z,C). 
Then the following are equivalent: 

(z,0 g WF(v^o)); 

(0, -^—) € WFseCe 4 ^) VW), t > t ; 
t — t 

(0, --^—) G WF sc (e i * < t . 

to — t 

This is a refinement of certain propagation results of which were in 
terms of quadratic scattering wavefront set WF qsc . In particular, one of the 
main results of JH] stated that if 8 = ir(Sb(z,C)) denotes the backwards 
limit of the geodesic through z, £ then for t > 0, 

(13) (0, -6/2t) i WF qsc Vo => C) i WF^(t). 

This clearly follows from Theorem [5] and the following wavefront set com- 
putations: 

Proposition 6. 

(9,0 G WF qsc u <=>(6,£- aO/2) G WF qsc e~ I£ff2 / 2 «. 

Proposition 7. // (0, 0) ^ WF qsc tt i/ien ^ WF sc u for any finite £. 

The proofs of Propositions EH3 which are not difficult, are best carried 
out using the more sophisticated definitions of wavefront set, involving as- 
sociated calculi of pseudo differential operators. 

8. Examples 

We conclude with a pair of examples which exhibit some extremes of 
behavior which solutions of Q may exhibit. 

First, we consider the Schrodinger equation in flat M. n with V = 0. Let 

Mz) = (-2™r 1/2 e~ l2?/2 ; 

note that this is just the one-dimensional fundamental solution, evaluated 
at time t = — 1, and extended to be constant in the variables z%, ■ ■ ■ , z n . 
The quadratic scattering wavefront set of ipo satisfies the hypotheses of (|13|) 
over two points in S™~ , given by (±1, 0, . . . , 0), hence (|T3|) guarantees that 
the wavefront set of ip(t) is confined at most to {t = 1, ( = (±1, 0, . . . , 0)}. 
On the other hand, we can compute exactly: z) = S(zi). Hence two 
points in the qsc wavefront set are able to produce an entire hyperplane of 
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singularities. By contrast, Theorem 03 gives a more satisfactory picture of 
the propagation phenomenon: we have 

(where z' = (Z2, ■ ■ ■ , z n )). Since the Fourier transform of this function is just 



i(n-l)/2^ 

we have 

WF sc (e- 2 /Vo) = {((±1, 0, . . . , 0), (0, C')) : C G K n_1 }- 
This is pair of hyperplanes over two points in S^" 1 which maps diffeomor- 
phically to WF^(l,z) according to the sojourn relation. 

Another extreme case for the results discussed above occurs if we work 
on M 1 with V = 0, and take 

M^)=e-' lz2/2 Ai(z). 

The quadratic-scattering wavefront set of tpo is just (—1,1/2) (Intuitively 
speaking, the factor of Ai kills the wavefront set at z = +1 owing to its 
exponential decay in that direction, but has no effect on the wavefront set 
of the e~ %z I 2 factor at z = — 1 owing to its slower oscillation.) Hence (|13|) 
permits wavefront set only at t = 1, in the direction £ = +1. On the other 
hand, we may compute exactly to find that 

il>(z, 1) = (-2^) 1 / 2 e i (^ + ^) E C°°(M); 

hence no wavefront set appears after all. This situation is accounted for 
in the results of ^S] by the use of a wavefront set that is uniform in time, 
and in terms of which this solution is singular everywhere on R at t = 1. In 
terms of Theorem^ which is about the wavefront set of ip(-,t) for each fixed 
t, the explanation is as follows: both e ir ^i/jq and ip[z, 1) have scattering 
wavefront set only at the corner S° x S°, alluded to in the remark following 
Definition [21 As there are no points (z,C) in WF sc (e* r ' 2 ipo) with £ finite, 
Theorem |S] asserts that ip(l,z) is smooth. 
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